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AFOSR  GRANT  #F4920-97-l-0497,  P00002 

(150-6756) 

“SHOCK  PROPAGATION  AND  SUPERSONIC  DRAG  IN 
LOW  TEMPERATURE  PLASMAS” 


SUPPORT  FOR: 


AFOSR-SPONSORED  WORKSHOP  ENTITLED: 
“Understanding  and  Control  of  Ionized  High-Speed  Flows” 
Princeton  University 
Princeton,  NJ  08544 
February  26-27, 1998 


ABSTRACT 

This  grant  supported  the  costs  of  an  Air  Force-sponsored  workshop,  “Understanding  and  Control 
of  Ionized  High-Speed  Flows,”  which  was  conducted  at  Princeton  University,  February  26-27, 
1998.  In  recent  years,  there  has  been  a  great  deal  of  interest  in  the  formation  mechanisms  and 
properties  of  air  plasmas.  This  interest  is,  in  part,  motivated  by  experiments  conducted  in  Russia 
and  in  the  United  States  which  indicate  shock  propagation  in  weakly  ionized  air  plasmas  is  at  a 
higher  velocity  than  would  be  predicted  by  presently  understood  models.  If  this  is,  indeed,  the 
case,  such  plasmas  could  be  used  for  supersonic/hypersonic  drag  reduction.  In  addition, 
atmospheric  plasmas  could  influence  flow  control  devices,  electromagnetic  attenuation,  and 
hypersonic  propulsion  systems.  As  a  consequence,  the  formation  of  such  plasmas  in  atmospheric 
pressure  environments,  and  the  study  of  the  properties  of  these  plasmas,  are  of  significant 
national  interest. 

DISCUSSION 

The  purpose  of  the  workshop  was  to  highlight  the  basic  science  issues  involved  with  controlling 
high-speed  flows  using  weakly-ionized  plasmas.  The  specific  objectives  of  the  two-day  meeting 
were  to  review  recent  experimental/theoretical  research  in  this  area,  and  to  exchange  ideas/ 
understanding  on  the  basic  flow  physics  of  the  problem,  and  to  discuss  fixture  research  directions. 
Princeton  University  was  a  natural  location  for  this  workshop  since  work  is  underway  here  to 
examine  both  drag  characteristics  and  electron  discharge  physics  associated  with  weakly  ionized 
plasmas. 

Sixty  individuals,  each  of  whom  is  considered  a  leading  authority  on  this  issue,  attended  the  two- 
day  workshop  (a  fiill  listing  of  the  participants  is  located  in  Appendix  A).  The  first  day  of  the 
workshop  consisted  of  four  sessions  on  an  agenda  of  twelve  speakers  and  covered  the  following 
topics: 


Session  I:  “Introduction  to  the  Problem”  (presentations  by  B.  Ganguly  of  AFRL  and  A. 
Auslender  of  NASA  LaRc). 

Session  H:  “Aerodynamics  of  Ionized  Flow”  (presentations  by  R.  Miles  of  Princeton 
University,  J.  Shang  of  AFRL,  N.  Malmuth  and  V.  Soloviev  of  Rockwell  Science,  and  V. 
Subramanian  of  The  Ohio  State  University). 

Session  HI:  “Theory  and  Computation  of  Ionized  Flows”  (presentations  by  G.  Kamiadakis 
of  Brown  University,  K,  Powell  of  the  University  of  Michigan,  and  S.  Nazarenko  of  the 
University  of  Arizona). 

Session  IV:  “Plasma  Generation  and  Maintenance”  (Presentations  by  S.  Kuo  of  Polytechnic 
University,  J.  Kline  of  Research  Support  Instruments,  and  W.  Rich  of  The  Ohio  State 
University). 

The  second  day  of  the  workshop,  the  participants  were  divided  into  four  separate  working  groups 
for  further  discussion,  followed  by  debriefings  and  laboratory  tours  at  Princeton. 

Appendix  A  is  a  complete  record  of  the  two-day  workshop  and  includes  the  final  agenda  for  the 
meeting,  the  complete  list  of  participants  and  their  contact  information,  as  well  as  copies  of  each 
of  the  presentations  for  all  of  the  four  topic  areas  listed  above. 
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FINAL  AGENDA 

AFOSR-SPONSORED  WORKSHOP: 


•’UNDERSTANDING  AND  CONTROL  OF  IONIZED 
HIGH-SPEED  FLOWS” 


February  26-27, 1998 

MEETING  LOCATION:  Convocation  Room  (C217)  Engineering  Quadrangle 

Olden  Street,  Princeton,  NJ  08544 


THURSDAY,  FEBRUARY  26, 1998: 


7:30  AM 
8:00  AM 


8:30  AM 

9:30  AM 
10:00  AM 
10:30  AM 


11:30  AM 
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1:30  PM 
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3:30  PM 
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DISCUSSION 
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Examine  Electromagnetic  Manipulation  Of  Ionized 
Flowfields  To  Improve  Hypersonic  Vehicle 
Performance 
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FINAL  AGENDA 

AFOSR-SPONSORED  WORKSHOP: 

"UNDERSTANDING  AND  CONTROL  OF  IONIZED 
HIGH-SPEED  FLOWS" 

February  26-27, 1998 

MEETING  LOCATION:  Convocation  Room  (C217)  Engineering  Quadrangle 

Olden  Street,  Princeton,  NJ  08544 


#  THURSDA Y,  FEBRUARY 26, 1998: 


7:30  AM 

SIGN-IN/COFFEE 

8:00  AM 

nPFNTMG  RFMARKS:  G.L.  Brown— Chairman,  MAE  Dept.,  Princeton 

R.  Miles-Princeton 

S.  Walker— AFOSR 

# 
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SESSION  I:  TNTRODUCTinN  TO  THE  PROBLEM 

B.  Ganguly- AFRL,  A.  Auslender-NASA  LaRC 
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DISCUSSION 
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10:30  AM 

SESSION  II:  AERODYNAMrCS  OF  IONIZED  FLQK 

R.  Miles-Princeton,  J.  Shang— AFRL, 
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N.  MalmuthA^.  Soloviev-Rockwell  Science 

V.  Subramanian— OSU 

11:30  AM 

DISCUSSION 

12:00  PM 

LUNCH 

1:30  PM 

CF<;«;Tn>J  tit-  THF.OR  Y  AND  COMPUTA  TION  OF  IONIZED  FLOWS 
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G.  Kamiadakis— Brown  U.,  K.  Powell— UM,  S.  Nazarenko— UA 

2:30  PM 

DISCUSSION 

3:00  PM 

BREAK 

3:30  PM 

QFSSTON  TV-  PT.ASMA  GENERATION  AND  MAINTENANCE. 

S.  Kuo-Polytechnic  Univ.,  J.  Brandenburg-RSI,  W.  Rich-OSU 

4:30  PM 

DISCUSSION 
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5:00  PM 

FIRST  DAY  WRAP-UP 

6:00  PM 

RECEPTION,  Prospect  House 

FRIDAY,  FEBRUARY 27, 1998: 
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8:00  AM 
8:15  AM 


10:30  AM 
11:00  AM 
12:00  PM 


SIGN-IN/COFFEE 
OPENING  REMARKS 
WORKING  GROUPS 
BREAK 

GROUP  DEBRIEFINGS 
LUNCH/LABORATORY  TOUR 
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Shock  Wave  Propagation  in  a  Weakly  Ionized 
Low-pressure  Discharge 


A.  1.  Klimov  et  al.,  Sov.  Tech.  Phys.  Lett,  8,  192  (1982) 


Shock  Stahdoff  Comparison  to  Predictions 

(0.75  in.  diam,  30  torr) 
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VS.  GAS  PRESSURE  AND  CURRENT  DENSITY 
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KINETIC  TO  ELECTROSTATIC  ENERGY  TRANSFER 


Weakly  Ionized  Gas  Research 

Russian  Ballistic  Tests  (G.  Mishin  &  A.  Kilmov,  1978) 


Schematic  of  Experiment 


Local  Valocity  in  Discharge  as  Function  of 
.  Discharge  Current 
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OUR  RLANNiD  RESEARCH  FOR  FY98-99 
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Long  range  cooperative  interactions  may  be  important 


A  Comparison  of  the  Shock  Thickness  With  a  0.8  mm 
Diameter  Wire  ,  30  Torr  Argon  ,  Pulser  Voltage  9kV 
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Argon  30  Torr,  10kV  2|iF,  42.2cm 
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El@ctrical  Sehematie 
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Deflectloiii  Signats  Downstream  of  Discharge 

Pulse  Energy  81  Joule 
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AFOSR  /  High-Speed  Workshop 


NASA  Langley  Research  Center 


AFOSR/ High-Speed  Workshop 
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Absorption  Ability  of  the  Microwave  Streamer  Discharge  j 
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The  Experimental  Data  about  a  Shock  Wave  Attenuation 

in  a  Microwave  Discharge 
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The  Acoustic  Interferometn*  Study  of  the 
Microwave  Discharge  Region 
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NASA  Langley  Research  Center 


Plasma  Drag  Reduction  Physics 


Development  of  related  mechanism(s)  and  analysis 
of  the  relevant  “fluid”  equation  set  (OSU) 
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Shock  dispersion  experiments  and  plasma  generation 

(future  work) 


“Anomalous”  Effects 
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B.  Supersonic  flow  in  weakly 
ionized  gas  around  a  spherical 
model  (ballistic  range 
experiments) 


Shock  Front  Velocity  in  Cylindrical  Glow  Discharge  [Phys.  Lett.,A230,  218-222  (1997)1 
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Ballistic  Range  Experiments 
Planar  Shock  Wave,  Positive  Column 
Planar  Shock  Wave,  Cathode  Region 


Relevant  Collision  Processes 
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Most  electrons  trapped  in  relaxation  zone 
Electrons  heated  by  inverse  bremstrahlung. 
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C.  Onset  of  turbulence  by  angular  momentum  transfer  in  energy 
pooling  reactions? 
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Precursor"  Microwave  Plasma  at  Mach  6 


Precursor"  Microwave  Plasma  at  Mach  6 


Microwave  Breakdown  in  Air,  O,  and  N 


p  (mm  Hg) 

Cw  breakdown  fields  for  air  and  nitrogen  at  a  frequency  of  994  Mc/sec. 
(Characteristic  diffusion  length  A  =  2.65  cm) 

A.  D.  MacDonald  et.  ai.  Phys.  Rev.  130, 1841  (1963) 


Microwave  Breakdown  in  Air,  O,  and  N 


(Characteristic  diffusion  iength  A  =  0.22  cm) 
A.  D.  MacDonaid  et.  ai.  Phys.  Rev.m  1841  (1903) 


PLASMA  DRAG  REDUCTION  MODEL 


nuclear  source 


PULSED,  DIGITAL  SCHLIEREN 


LASER/MICROWAVE  INTERFEROMETRY 


PULSED,  DIGITAL  SCHLIEREN 
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Program  Overview 

Shock  Wave  Dynamics  in  Weakly  Ionized  Plasmas 

at 

or  Modem  Fluid  Physics,  Florida  A&M  University  (FAMU),  Tallahassee,  FL 
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duction  of  operational  aircraft. 


Recent  FAMU  Achievements 
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Modifications  of  our  current  Supersonic  Shocked  Plasma  Facility  to  afford  a  wider 
variety  of  manipulations  in  the  shock  wave-plasma  interaction,  including  the 


troduction  on  operational  airframe  platforms. 


AERODYNAMIC  APPLICATIONS 
OF  WEAKLY  IONIZED  PLASMAS 

R.  Miles,  S.  Macheret,  S.H.  Lam,  R  Efthimion,  L  Martinelli 
Princeton  University 


Drag  Reduction 
Plasma  Ramparts 

Flow  Control 

By  Localized  Energy  Addition 
By  J  X  B  forces 

Control  and  Intiation  of  Chemical  Processes 

Power  Extraction 

MHD  in  an  engine  inlet  (AJAX) 

MHD  in  ionized  flow 

Hypersonic  Wind  Tunnel  Drivers 

By  Energy  Addition  through  microwaves 
By  MHD  acceleration 

Suppression  of  shock  strength  and  sonic  boom 


SUSTAINING  A  PLASMA  IN  HIGH  SPEED  FLOW 


Important  for  Drag  reduction.  Plasma  Ramparts,  Flow  Control, 
Shock  Suppression,  MHD,  and  Combustion  Initiation  and  Control 


Discharge  may  be  blown  out 


Thermal  (equilibrium)  plasma  will  need  to  be  reinitiated  - 
a  thermal  plasma  cannot  propagate  at  supersonic  speed 


Nonequilibrium  plasma  will  require  high  electric  field  and  high  power 


A  filamentary  plasma  may  serve  the  purpose  and  will  require  less 
power  (a  small  portion  of  the  volume  is  ionized) 


Filaments  may  be  guided  by  laser  or  electron  beams 


A  high  power  microwave  may  be  the  preferred  source. 


SUSTAINING  MICROWAVE  DRIVEN 
PLASMAS  IN  ATMOSPHERIC  PRESSURE  AIR 


Thermal  Regime 

Requires  >  4000K 

Needs  relatively  low  microwave  power  (kilowatts/cm“) 
the  power  offsets  the  thermal  losses 

Long  lifetime  (cooling  time) 


Nonequilibrium  Regime 

Requires  high  power  microwaves  to  overcome  the  electron 
recombination  and  attachment  rates 

Diffuse  plasma  (low  pressure) 

Low  temperature  neutral  species 


Filamentary  discharge  (high  pressure) 

High  temperature  localized  streamers 


FILAMENTATION  OF  ATMOSPHERIC 
MICROWAVE  DRIVEN  PLASMA 
Grachev,  Esakov,  Mishin  and  Khodataev 

Moscow  RadioTechnical  Institute 


LASER  GUIDED  ARC  FORMATION 


ELECTRODE  : 

-  •  ,  ■  Tv  - 

i  ^  ;^&ECTRODE 
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drag  reduction 

Expciinicntcil  Obscrvcitions 


Shocks  Accelerate  Upon  Entry  Into  Plasmas  From  Neutral  Gas 

Dramatic  Changes  In  Structure; 

Shock  Front  Broadens 

Amplitude  Decreases 

Shock  Typically  Splits  Into  Two  or  Three 

Effect  Persists  for  Quite  a  Long  Time  (up  to  1  - 10  msec) 


TrtXS’L  (Longitudinal  And 

nsverse),In  Various  Gases  (Air.  Argon,  Hydrogen.  Helium  Etc  1 
Using  Various  Diagnostics.  nuium,  ntc.) 


Ballistic  Range  Experiments  Seem  to  Imply  a  Chance  in  the 
Effective  Speed  of  Sound,  and  Show  Decrease  in  Su^c  Drag. 


DE 


Planar  C 


Planar 


EXPERIMENTAL  CONCERNS 


Temperature 

In  all  the  experiments,  the  temperature  of  the 
neutral  gas  has  not  been  directly  measured. 


Temperature  Profile 

A  transverse  temperature  profile  or  a  curved 
thermal  boundary  will  lead  to  complex  shock 
structure.  This  may  occur  as  the  shock  passes 
through  the  hot  cathode  or  anode  fall  regions 
of  the  discharge.  (Well  away  from  the  centerline 
of  the  experiment) 

Models  indicate  the  apparent  shock  splitting 
and  the  density  rise  echo  experimental 
observations 


Nonuniform  Temperature  Field 

Shock  splitting  has  been  seen  in  turbulent 
mixtures  of  gases.  If  the  plasma  has  hot  and  cold 
regions,  a  similar  shock  break-up  may  occur. 


CONCLUSIONS  ON  DRAG  REDUCTION 


All  the  experimental  phenomena  we  examined  are 
consistent  witli  thermal  gradient  effects 


Ballistic  range  data  needs  to  be  studied  in  more  detail 


Accurate  temperature  and  density  measurements  need 
to  be  made 


There  may  be  some  effects  of  thermal  gradients,  unsteady 
phenomena,  vorticity,  and  random  or  flickering  heat  addition. 


PRINCETON  PROGRAM 


Development  of  new  diagnostics  for  temperature  and 
density  imaging 


New  World  Vistas 

Sustaining  microwave  driven  plasmas  in  a  Mach  3  flow 
Plasma  effects  on  the  bow  shock  in  a  Mach  3  flow. 
Study  of  basic  mechanisms  of  acoustic  wave 
interactions  and  shock  dynamics  in  stationary  plasmas 


Plasma  Ramparts 

Experimental  studies  of  the  control  of  plasma 
filamentation 

Kinetic  model  development  for  pulsed  ionized 
air  plasmas 

Modeling  nonequilibrium  molecular  processes 
High  power  microwave  driven  discharge  dynamics 
Stabilizing  and  extending  microwave  driven  plasmas 
with  magnetic  fields 


Mariah  II  (Radiatively  Driven  Wind  Tunnel  augmented  by  MHD) 
Sustaining  conductivity  in  atmospheric  pressure  air 
with  electron  beams 

Acceleration  and  deceleration  of  supersonic  flows  in 
MHD  channels 

Energy  couphng  into  high  pressure  air 
MHD  acceleration  with  filamentary  discharges 


PLASMA  TEMPERATURE  MEASUREMENT 
BY  FILTERED  RAYLEIGH  SCATTERING 


Rayleigh  Scattering 

Total  Scattering  is  Proportional  to  Density 

Linewidth  is  proportional  to  Jt 
Signal  is  weak 

Background  scattering  is  at  nearly  the  same 
frequency 

Filtered  Rayleigh  Scattering 

Eliminates  Background  Scattering 
Permits  Measurement  of  Linewidth 


t  Laser 


A 

b 

Rayleigh  scattering  +  background 

c 

filter  profile 

d 

imaged  signal  (no  background) 

FILTERED  RAYLEIGH  SCATTERING 
IMPLEMENTATION 


Source 


Filter 


Camera 


Tunable  Ultraviolet,  Injection  Locked  Laser 
Ti: sapphire  with  “Ramp  and  Lock”  Technology 
(Schwartz  Electro-Optics) 


Atomic  Mercury  Vapor 

2”  Diameter  Cell  Held  at  40°C  with  a 
hot  water  bath  (vapor  pressure  1  Torr) 


UV  Sensitive,  Intensified  CCD  (Princeton 
Instruments) 


SEQUENTIAL  RAYLEIGH  IMAGES  OF  MACH  2.5 
SHOCKWAVE  /  BOUNDARY  LAYER  INTERACl'lON 

AT  A  14°  WEDGE 


Framing  rale  is  500,000  images  per  second 


Mariah  II 


No  wind  tunnel  can  at  present  be  run  to  accurately 
simulate  flight  conditiona  at  higher  than  Mach  8 

240K,  (120,000  ft)  requires  2800K,  8  atm  in  plenum 
NO  contamination 

Low  static  T  or  low  static  P  (enthalpy 
too  low  or  entropy  too  high) 

Short  time  operation 

Radiatively  Driven  concept  adds  energy  downstream  of 
the  throat  in  the  supersonic  section 

Static  temperature  stays  low 
Low  temperature  but  high  pressure  in  the  plenum 
Concept  is  driven  by  lasers  or  electron  beams 
Demonstrated  with  lOKW  laser  in  December  at 
Wright  Labs 

MHD  section  is  to  be  added  to  further  accelerate  the 
flow  after  the  Radiative  section 

Reduces  the  front  end  pressure  requirement 
Extends  the  envelop  of  the  tunnel 
Requires  electron  beam  sustained 
conductivity 

Must  operate  at  close  to  room  temperature 
and  at  pressures  on  the  order  of  .1  atm. 


Energy  Addition  by  Lasers 


RDHWT  Review  April  23  -  24,^1997 


Figure  1.  Diagram  of  10  kw  laser  energy  addition  experiment 
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Thermodynamic  limit  of  RDHWT: 

Mi=T^As 

Ta  can’t  be  extremely  high  (erosion,  chemistry); 
minimum  ^  is  defined  by  incompressibility; 
maximum  s  is  defined  by  test  section  requirements. 

THERMODYNAMIC  ADVANTAGE  OF  MHD 

Total  enthalpy  augmentation: 

Ah=-^TAs 

‘  K-l 

where  K^EI(uB)=  (Energy  added)/(Push  work) 

The  best  thermodynamically: 

when  Aht  can  be  large  even  for  small  T  As 

(non-entropy-generating  energy  addition). 

However,  when  A->1,  j  =  (j  {E-uB)  ->0, 
and  to  add 

Ah,  =  GK{K-l)uB^L/p, 

very  large  conductivity,  a,  and/or  length,  i,  would  be 
required 


AJAX 

Proposed  by  Russians  in  1990 
(Vladimir  Freistadt  - 

State  Hypersonic  Systems  Scientific  Research  Enterprise 

St.  Petersburg) 


MHD  converter  in  the  inlet  connected  to  MHD  accelerator  in  the 
nozzle  to  by-pass  free  stream  kinetic  energy  and  reduce  the  entropy 
rise  in  the  inlet. 

Reduces  the  flow  to  subsonic  in  the  combustor 
MHD  conditions  almost  the  same  as  Mariah  II 
Requires  some  method  to  sustain  conductivity 


Reduction  in  vehicle  drag  by  cold  plasma  injection  at  the  nose 


Increase  of  combustion  volume  and  efficiency  within  the  engine  by 
plasma  injection  or  injection  of  materials  ahead  of  the  fuel  injectors 


Reforming  of  kerosene  and  water  into  hydrogen  and  CO  for  fuel. 
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Design  new  experiments 

Compare  evoiving  theoreticai  modeis  with  previous  datasets 


Supercruise  Fighter 

Plasma  Aerodynamic  Enhancements 
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-  hypersonic  applications 
P4  study  useful  since 
experimental  dataset 
available  to  test  models 


Problems 


Problem  3a  -  (F-15) 


1-D  shock  propagating  into  a  plasma 


Perfect  Gas  Shock-Layer  Problem 


S  s  viscous  shock  layer  thickness 
-  ]']  s  shock  strength 
Y  s  specific  heat  ratio 


Scaling  of  Shock  Attenuation  with 
Thickness  in  a  Piasma 


Drag  Reduction  Mechanisms 
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Flow  Configurations 


Aerospike  Effect  on  Pressures 


Peak  attenuates,  with  increasing  spike  length  reducing  drag 

Needs  to  be  traded  against  high  heat  transfer  near  peak 

PMma!  heating  can  affect  reattachment  and  therefore  pressure  peak 


Estimate  of  Drag  Reduction 
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Cumulative  Pressure  and  Skin  Friction  Drag  with  Plasma  Heating 


(MM) 


Comparison  of  Predicted  Drags  With  Russian 
Experiments 
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Forebody  drag  of  a  cone-cylinder  with  counterflow  jet  in 
supersonic  flow _ _ _ _ 


Science  Center 


Counterflow  jet,  =  2.0,  T^jT^  =  10,  7  =  1.4. 


Forebody  drag  of  a  cone-cylinder  with  counterflow  jet  in 
supersonic  flow _ ■  _ _ 
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vorticity,  multiple  secondary  shocks  and  instabilities 


Forebody  drag  of  a  cone-cylinder  with  counterflow  jet  in 
supersonic  flow 


Science  Center 
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Figure  7:  Forebody  drag  coefficient  as  a  function  of  PaIPoo  for  different  values  of  TaIT^o  as 
follows:  crosses:!,  triangles:3,  diamonds:6,  asterisks:10.  The  forebody  drag  coefficients  of  the 
sharp  cone  and  the  cone  frustum  are  shown  as  lines  for  reference.  Mo©  =  2,  7  =  1.4 
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Figure  8:  The  minimum  value  of  the  forebody  drag  coefficient  as  a  function  of  Ta/T^. 
■Woo  =  2,  7  =  1.4 


Propagation  of  1-D  Shock  Waves  Through  Weakly  Ionized 

Plasmas 

V.R.Soloviev,  V.M.Krivtsov,  A.M.Konchakov 
Moscow  Institute  of  Physics  and  Technology,  Russia 

N.Malmuth 

Rockwell  Science  Center,  CA,  USA 

Main  topics  of  the  model 

•  Hugoniot  relations  on  a  shock  wave  (SW)  in  a  weakly  ionized  plasma 

•  Estimations  of  the  energy  stored  in  different  excited  states  for  glow  discharge 
plasma 

•  The  existence  of  high  situated  metastable  states  in  Ar  and  N2 ,  which  serve  as 
a  ‘new  ground  state’  for  higher  states  excitation  by  electron  impact 

•  Hi^ly  excited  atomic  and/or  molecular  states  store  sufficient  energy  to 
produce  SW  acceleration  and  attenuation 

•  Inclusion  of  highly  excited  atom  energy  release  behind  the  SW  front 
dramatically  improves,  agreement  of  tlieory  with  experiment  for  1-D  shock 
penetrating  into  a  plasma  and  hypersonic  sphere  entering  a  plasma 


Hugoniot  relations  on  a  shock  wave  (SW)  in  a  weakly  ionized 

plasma 


A  set  of  gas  mass,  momentum  and  energy  conservation  equations  in  SW 

coordinates 

dt 


- —+ — (oit(u  +  D)  +  p)  =  — 

dt  ^  d4 
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p,  p,  h  are  density,  pressure  and  enthalpy  of  the  gas;  v  is  its  velocity; 

^  ^  are  thermal  conductivity  and  dynamic  viscosity  coefficients; 

Cp,  j'are  specific  heat  at  a  constant  pressure  and  the  ratio  of  specific  heats; 
D  is  the  SW  velocity; 

The y^and  Or  terms  describe  the  plasma  influence  on  the  gas. 


Hugoniot  relations 


p,u^=p^D, 
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To  produce  a  plasma  effect  on  SW  propagation  it  is  necessary  Sui,  Su2,  or  eu  to 
be  of  the  order  of  1. 


The  typical  cold  glow  discharge  plasma  parameters  in 
SW  -  discharge  interaction  experiments 


gas  pressure  (gas  density) 
gas  temperature 
current  density 
electric  field  strength 
E/N  value 

discharge  specific  power 
electron  temperature 
electron  density 

metastable  atom  (or  molecule)  density 


p  «  3  -  30  Torr  ( N  «  10^*^  -  10^*  cm'^ ) 
T«300K 
j  «  1  - 100  inA/cm^ 

E  « 100  V/cm 
(2  -10)  xlO*’^  Vxcm- 
W«0.1  - 10  W/cm^ 

Te «  1  -  3  eV 
lie®  lO'®  - 10*"  cm'^ 

10**  - 10*^  cm'^ 


The  energy  stored  in  different  excited  states 

•  freestream  gas  entlialpy  ho=NcpT  ~  10*^  eV/cm^ 

•  ionization  niEi  ~  10*^  eV/cm^,  Sh=  iiiEi/NCpT  ~  10“^ 

•  metastable  atoms  (or  molecules)  NmEm-  lO*"*  eV/cm^,  8h=  NmEm  /NcpT~10'^ 

•  specific  discharge  energy  supply 

during  the  time  of  unit  length  S  W  pass  Eh=  Wx  1  cm/D  NcpT~l  0'^  - 1  O'* 

•  the  vibrational  energy  is  shown  to  be  imsuitable  to  explain  the  SW  anomalious 
propagation  effects  because  of  large  vibrational-translational  relaxation  time 
TvT~lfi 


The  scheme  of  Ar  and  N2  electronic  excited  states 


A\]5.7^\) 


N2(15.58ey) 


4p(B2) 
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Ar(4s)  and  N2(A^2)  are  the  metastable  states 


The  At  metastable  state  quenching  and  building  up 


Process 

Ar(4s)  +  2Ar  ->  Ar*2  +  Ar 
Ar(4s)  +  e->-  Ar  +  e+1  L8ev 
Ar(4s)  +  N2  ->  At  +  NsCC,  B, ...) 
Ar(4s) +Ar(4s) ->  Ar"^  +Ar  +  e 
Ar(4p)  +  Ar  ->Ar(4s)  +  Ar 
Ar(4p)  ->Ar(4s)  +  hv 


Rate  constant  Quenching  time  for 

p=10  Ton,  ne=10^^cm*^ 
and  1%  N2  diluent  in  Ar 


10'^^  cmVs 

0.8  ms 

10'^  cmVs 

1ms 

3.6xlO'”cmVs 

0.01ms 

10'^  cmVs 

1  msx(  1 0  ’■/n(4s)[cm' 

~  10'”-  10*’"cmVs  0.3  -  3mcs 

3.2x10^  s'’ 

0.03mcs 

The  N2(A^X)  metastable  state  quenching  and  building  up 


N2(A^Z)  +  Ar  ^  N2(X)  +  Ar 
NzCA^S)  +  N2(X) ->  2  N2(X) 

N2(A^Z)  +  O2  ->  N2(X)  +  O2 

N2(A^S)  +  N2(X,v)  N2(B^n)  +  N2(X) 

N2(A^2)  +  N2(A^S)  ->  N2(B,C,...)  +  N2(X) 
NiCB^n)  +  N2(X)  NzCA^Z)  +  N2(X) 
N2(B^I1) ->N2(A^Z)  +  hv 
N2(C^n)  ->N2(A,B)  +  hv 


<4x  1 0‘*'*cm^/s  >  0. 1ms 

<4.5xlO’'^cmVs  >100ms 

3.6xlO'’^cmVs  4mcs 

(for  air  composition) 

~  10"^°cm^/s 

1 0'^cm^/s  1  msx(  1 0*^/n(A)[cm'^]) 

2x10*^ ’em  Vs  O.lSmcs 

1.5x10^  s*’  6.5mcs 

2.7x10*’  s'’  0.037mcs 


The  kinetic  model  of  energy  storing  in  highly  excited  states 


The  system  of  balance  equations  for  excited  states 


^  +-i-)  -  k,n,N,  +  N\k^N  +  A.ff) 

m  q 

^N* 

-^  =  Kn.N„-N\kN  +  A^e) 


E-E. 


exp( — - for  n>n*  kjie»k^ — 


ErE*^0.5-UY 

%  e  e  e 


Ry  =13.6eV  is  the  Rydberg  constant 


W=jE  -  is  the  discharge  specific  power; 

Ti~l  -  is  the  discharge  fractional  power  input  into  excitation; 

Em,  En-  are  the  energies  of  the  different  states  excitation  from  the  ground  state; 

-  is  the  metastable  atom  Penning  process  rate  constant; 

Tq  -  is  the  metastable  state  quenching  time  in  other  processes; 

ke  -  is  the  electron  excitation  rate  constant  from  metastable  to  N*  states; 

kq  -  is  the  gas  quenching  rate  constant  for  N*  states; 

Ar-  is  the  radiative  decay  rate  for  N*  states; 

0  -  is  the  probability  of  radiation  trapping; 

AEi  -  is  the  decrease  of  ionization  potential  due  to  excited  atom  electron  -  ion  interaction; 


For  the  tube  radius  R 
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e  = 


For  Ar(4p)  Ar(4s)  transitions  0f«3xlO’^  (10^^/Nm)^  ^  followed  by 

AR0«kqN 


The  steady  state  values  of  Nm  and  N* : 


I— ^  «  25  X 1 0"  y]w[W /  cm^],  cm~^ 


N  «  3.5x1 - -,cm^ 
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k,N 


/?[7brr] 


The  equilibrium  highly  excited  hydrogen-like  states 
population 


n 
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-  the  upper  limit  of  principle  quantmn  number  for 
hydrogen-like  states; 
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g*  ~30  is  tlie  degeneracy  of  N*  state 


The  energy  stored  in  highly  excited  states 
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The  energy  released  to  heat  behind  the  shock  wave 

N  N 

0  and  1  indexes  correspond  to  the  gas  and  plasma  parameters  before  and  behind 
the  shock 


The  relaxation  length  behind  the  SW  front 
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Inclusion  of  Plasma  Effects  Dramatically 
mproves  Agreement  of  Theory  with  Experiment 
for  1-D  Shock  Penetrating  Into  a  Plasma 
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Inclusion  of  Plasma  Effects  Dramatically 
mproves  Agreement  of  Theory  with  Experiment 
for  Hypersonic  Sphere  Entering  a  Plasma 
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Conclusions 


•  The  mechanism  of  shock  wave  (SW)  acceleration  and 
attenuation  due  to  interaction  with  a  weakly  ionized  plasma 
really  exists 

•  Preliminaiy  calculations  of  1-D  SW  speed  in  plasma  are  in  a 
good  agreement  with  experimental  results  when  the  highly 
excited  states  energy  release  behind  the  SW  front  is  taken 
into  account 

•  SW  acceleration  is  accompanied  by  decrease  of  density  and 
velocity  jump  with  a  small  increase  of  static  pressure  jump 
through  the  SW  front 

•  The  foregoing  results  in  a  decrease  of  dynamic  pressure 
jump  through  the  SW  front 

•  The  creation  of  a  discharge  plasma  IN  FRONT  OF  the  shock 
wave  courses  the  SW  attenuation  and  drag  reduction  for 
supersonic  and  hypersonic  vehicles 

•  More  detailed  kinetic  analysis  is  necessary  to  obtain  more 
precise  results 


Non-Equilibrium  Thermodynamics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 
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Normalized  Density  Normalized  Density 
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—  shock  nearly  recovers  its  original  strength 
well  downstream  of  the  discharge. 
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but  cannot  explain  others  (such  as  shock 
recovery  distance). 
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Suppose  that  space  charge  layer  induces 
a  potential  in  the  neutral  gas,  because 

neutrals  are  polarizable: 


<  (2e\|//m) 


p 


I 


V 


Non-Equilibrium  Thennodynainics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 


o 

O 

U 


0) 

G 

Gu 

0) 


•  • 

o 

0) 

o 

JGl 

c/) 

QJ 

pfl 

"o 

•  ^m^ 

Vi 

G 

a 

o 


•  A 

O 

o 

e 

Co 

•p<k 

•r^ 

C3 


O 

^C) 

o 

|s 

o 

o 


o 

o 

V 

A 

N 

N 

.Co 

r\ 

c\ 

•k.4 

O 

T— 1 

^  K 

■■  '  ■  - 

II 

O 

N 

n 


particle  velocity 


Non-Equilibrium  Thennodynamics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 


•pM  ^ 

ng  O 

25  ^ 

S  ^ 
S  ^ 

O 

— 

^  o 
o  ^ 


(N 


(N 

(N 


E-s 

K 

(N 

I 

+ 

CO  ^ 


February  26, 1998 


nil! 


Energy 


Source  terms  in  continuity  equation 
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Flow  Control  via  MHD  and  EMHD 

George  Em  Karniadakis 

Division  of  Applied  Mathematics 
Center  for  Fluid  Mechanics,  Turbulence  and  Computation 

Brown  University 
Providence.  RI 02912 

http://www.  cfm.  brown.  edu/CR  UNCH 
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Understanding  and  Control  of  Ionized  High-Speed  Flows 
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High-Speed  Flow  Features 


Strong  shock  wave  effects/interactions 

-  thin  layers,  small  detachment  distance,  etropy  layer 


Strong  viscous  effectsAnteractions 


-  increased  drag,  heating  at  leading  edge 

Transitional  flow 

-  Rtc  ~  10*  vs.  10^  (subsonic) 

-  Cj  and  Ch  are  3x  laminar 

-  Rti  ~  lOOMe 

Nonequilibrium 

-  thermal,  chemical 

-  introduction  of  length  scale 

-  multi-rate  physics 

Ionization 


-  4000  K  to  6000  K:  NO  NO"^  +  e"  (mild) 


-  >  9000  K:  +  e- 

0  —*■  O’*"  -b  c“ 


-  Above  20kM  +  10  km/s 


-  External  E.  frequencies 

Multiple  boundary  layers 


-  hydrodynamics,  thermal,  species 

-  electric  conductivity,  etropy  layers 


Non-Equilibrium  Plasma  Boundary  Layer 

-  Effects  on  Te,  cr 

-  atmospheric-pressure  potassium-seeded  arpon  plasma 


Normalized  Temperature  Normalized  Conductivity 

Ref;  Brown  &  Mitchner  (1971) 

Even  in  plasma  which  is  strongly  collision  dominated,  significant  nonequilibrium 
may  be  present  in  the  boundary  layer  adjacent  to  a  cooled  solid  surface...” 


Note:  Joule  heating  due  to  E  could  further  increase  conductivity. 


Electrical  Conductivity 


Electric  Field  vs.  Ionization 


•  In  constant  E,  B  fields  tensor  conductivity 

'  0^11  =^(nc,me,I/eF,/°) 

Weakly-ionized 


me,  I/ei/,  f,  UJe) 

i  me,  I'eHi  (*^e) 


-  If  B  =  0  aj^  =  =z  a 


-  =»  dependence  of  a  on  a;  ^  ^  U»'\.  | 

Note:  /®  =  ,  Ex,  1/4^^,  T47,  m*,  m^f ) 


Lorentz  Force:  Numerical  Modeling 


•  Governing  Equations  for  B  and  E 

•  Boundary  Conditions  -  The  Motz  problem 
Example  1:  Electromagnetic  Tiles 


Diiichlet: 


;  Neumann: 


F±_  jP|| 


Example  2:  Alternate  stripes  of  electrodes/magnets  e 


Dirichlet  Neumann  D-N  Mixed 


Lorentz  Force 


•  MHD  Assumptions  {tc  >  a; 


-  Convection  current  neglected:  — —  <C  1 

nJ 


Force  due  to  net  charge  neglected:  -j-j.  -  — j  -C  1 


Then: 


Fl  —  J  X  B\  Eo=  external 
J  =  (7{E  +  uxB)  +  aEo  Okn'^ 


MHD-Control 


EMHD-Control 


Apply  B 


•  E  =  -V$ 


Apply  B  and  E 


•  J  cfEq 


V2$  =  V.(uxB)  • /applied  = 


(tEqBL 


aB^L 


‘induced 


where  R^rn  =  iicrUL  <  1  in  lab  scale  ^Ac^Uc;t  fi'e(As^ 


•  But  Ohm’s  law 


tiallv  ionized  gas: 


J  +  /?e J  xb  +  sb  X  {J  X  6)  =  a(E  +  ux  B)  +  Vpc 

1^  CTlg 

j  “Generalized  Ohm’s  Law” 


Discontinuous  Galerkin  Method  -  MskT olt 


dynamical  Direct  Numerical  Simulation  ^DNS) 

MHD,  Compressible  &  Incompressible  flows 

Full  3D  Configurations  -  Arbitrary  Geometric  Complexity 

Unstructured/Hybrids  Grids 


High-Order/Spectral  Accuracy 

Conservative  Formulation 

h-refinement  for  shocks/discontinuities 

ALE  Algorithm  for  Moving  Rotating  Subdomains 
METIS/MPI-based  Parallelization 


K4CH2./S 


H-P  Adaptive  Refinement:  Viscous  Supersonic  Flows 


Figure  1:  Discontinuous  Galerkin  Simulation:  Density  contours 


Figure  2:  Unstructured  grid  and  variable-order  per  element 


Jt  APi^e-oAc-vf  : 

•DKlS  -  <Neo.r 

•  L-f  l'<\  c<^u4Ufc  AtJtiA  iT^^eAi 


Hybrid-Element  Boundary  Layer  Resolution 


0.0  0.5  1.0  1.5 


PhD  Thesis:  T.Warburton 


fJDAjS  li'j>  /C,f>ne.A*<-w.'('- 

^  Ja  equ«:Ut>/;aA  lO.U 

<4 

^  Ui  ok  ^  ,  CjQt^^f  kL^ 


MHD/EMHD-Simulations 


1.  Laminar  &  Turbulent  Wakes 

-  vortex  street  suppression 

-  phase/frequency  control 

n.  Wall  Turbulence 

-  High  <7,  spanwise  forcing/MHD 

-  Low  <7,  steamwise  forcing 
0  -  Low  (7,  “normal”  forcing 

Method:  Parallel  DNS  -  The  Jsfen'T ocr  Code 


Experiments:  Conflicting  results 


f  FZR  TU  Dresden 


T.  Weier,  O.  Posdziech 


EMHD  control  of  the  flow  around  a 

circular  cylinder 


i  .  ■■  - . . 1  ^ - - - - - - - 

Dimensionless  numbers: 

-  Reynolds  number 

jy  _  JO  Bq  d 

-  Interaction  Parameter 

P'^lo 

Air  Force  Workshop 


Princeton  NJ 


Feb.  25,  1998 


Workshop  Princeton  NJ 


TJR>aocg>^T  vxy^fegx  4r4A- 


Cylinder,  Re  =  500 


Active  Control:  Lorentz  Body  Force 


»  Ph‘l'  ^  - ►  Electric  Field 

•  ^  Magnetic  Field 

»  4h  - ►  Lorentz  Force 


^  ^  OE 

4  <kro«.R. 


Nosenchuck  &  Brown,  1994 


Summary 


•  Explore  MHD,  EMHD  and  EHD 

•  Strong  electrical  conductivity  non-uniformity 

•  Selection  of  physical  model  is  application-dependent 

-  degree  of  ionization 

•  Transitional  Flow  -  No  arf  hoc  models 

•  Numerical  Model  (A/* e/cT cxr  ): 

-  High  order 

-  Hybrid  grids 


•  Benchmark  Experiments 
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Spark  Morphology 

Velocity  Contours 


Left  Boundary:  p=i.o;  d  =  i.o;  u  =  v  =  o 

Right  Boundary:  p  =  O.l;  d  =  "V";  =  v  =  0 
(top  and  bottom  B.C's  are  "reflecting  slip") 

Initially/  these  two  states 
meet  at  a  discont  inuity* 


DENSITY:  "Y" 


B>arocUn\c  \/of4\o4^ 


CO  =  VKH 


- ((^■v)u  =  -^v-ii-Vj y-vf> 


Ixirocl'ii^lc 


Po4en4;^l  Vorficif^  ConS«VA4ioi» 

2  J)  bcxT of rop  I  c  (p=p(s'>)  -P^U 


Potenfift^  Vorfjci'K  ^ 
eocK  -P-Cw-xti.  e?i2.ir»aK  ■ . 
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STMULATIOM  METHOD; 

(of  the  2D  Euler  equations) 

PHYSICAL  ASPECT  RATIO  (x:y)  =  100:1 

#  of  gridpoints  in  x  direction  =  400 

#  of  gridpoints  in  y  direction  =  20 

#  of  ghostpoints  =3  (on  each  side) 


1:  Given  tlie  fluid  parameters 
at  each,  gridpoint,  find  th.e 
eigenvectors  to  Roe's  matrices* 

2:  Use  a  Sth.  order  ’weigh.ted  ENO 
sclieme  to  find  t  lie  flux 
eigenvectors  at  Ixalf  gridpoints  • 


3:  Use  a  3rd  order  Runge  Kutta  routine 
to  propagate  tlie  lialf  gridpoint 
fluxes  and  find  tlie  nev  fluid 
parameters  at  tli  e  real  gridpoints  • 


4 :  Repeat • 
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RSI  Further  Investigation  of  Large  Volume  PI  A 

Research  (Persistent  lonizationin  Air)  Plasmas  at 

Support  Atmospheric  Pressure 


Materials  and  Methods  2.45  GHz  (1) 


RSI 

Research 

Support  Materials  and  Methods  2.45  GHz  (3) 
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Plasma  Torch<es  and  Their  Demonstration 


Spencer  P.  Kuo  and  Edward  Koretzky 


Department  of  Electrical  Engineering 
Polytechnic  University 
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I  t-4 
CJ  *3 
O)  CO 

"Z  a 


d 

g  tS  g 

Q3  O 

•s  g :: 

"si 

<ii 

d  Cu 

•  *  1-4 

J-4  t4-l 

^  d  O 


<D  (U 
CO  ^.3 

d  cd 

03  ^ 

S  z 
$  «2 
"O  CO 


o  ct 


O  "g  CO 
II4  QJ 

?n  H  X 

hO 

r-r  ^  CO 

U4  U  > — ^ 


03  CL 
■q3  cy 

y, 

a  o 

4-9  Ui 

CO  4-9 


jet  nozzle. 


Plcisma  Torch 
1  ATM 


Pressure: 

Density: 

Temperature: 

Dimensions: 

Volume: 

Power  Consumption: 

Peak  Voltage  and  Current: 


1.5  X  10^^  (deduced  from  mea¬ 
surements  based  on  two  inde¬ 
pendent  methods) 

1200"  A' 

cylindrical  shape  having  1  cm 
radius  and  7  cm  length 

20  cm^/torch 

600  W/torch 

2.8  kV  and  2  A 


Microwave  Plasma  Layers 

Pressure:  20  torr 

Density:  10^^  cm~^ 

Microwave  Power:  1  MW  peak,  1  fis  pulse,  40  rps 
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Plasma  Ramparts  Using 
Metastable  Molecules 

MURI  Program 

The  Ohio  State  University 
Princeton  University 

"Optically  Pumped  Nonequilibrium 

Plasmas" 

L  Adamovich,  V.  Subramaniam,  W.  Rich 

Presented  at 

AFOSR  Workshop 
"Understanding  and  Control  of 
Ionized  High-Speed  Flows" 

Princeton  University 
Feb.  26-27, 1998 

Work  funded  by  the  Direetor  of  Defense  Research  & 
Engineering  (DDR&E)  within  the  Air  Plasma  Ramparts 
MURI  Program  managed  by  AFOSR 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Program  Objectives: 

•  Create  large  volume  (0[m^])  free  air  plasma 

•  One  atmosphere 

•  Free  electron  densities  10^^  cm^  or  greater 

•  Gas  temperature  less  than  2,000  K 

•  Energy  efficient  method 

Approaches: 

A.  Thermal  Plasmas  (High  Temperature  Arcs,  Reentry 
Plasmas) 

Problems: 

I)  High  Temperature 


B.  NonThermal  Plasmas  (Glow  Discharges,  C 
Pumped  Plasmas) 


Problems: 

I)  Requires  thermodynamic  work  (electrical  work, 

laser  work) 

II)  Stability 

OSU/Princeton  Approach:  Laser  Pumped  Metastable 
State  Plasmas. 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


92  W  cw  broad-band  lasing  conditions 


Date:  10-26-1992 _ -i 


Nonequilibrium  Thermodynamics  Laboratories 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


The  Ohio  State  University 


Advantages: 

Stable,  large  volume 
High  pressure 
Low  gas  temperatures 

Problems: 

Electron  density 
Efficiency 
Non-air  species 


Nonequilibrium  Thermodynamics  Laboratories 


Three  Level  Laser 
Energy 


Pumping  Waste  Heat 


Upper  level  energy:  Eu 
Lower  level  energy:  El 


Ground  level  energy:  Eo-0 


Quantum  Efficiency 
_  Eu  -  El 

^  qudntum  g 

L 

-40%  forC02/N2 
-95%  forCO/N2 


Scoville,  1954  showed  equivalence 
with  reversible  engine  efficiency 


The  Ohio  Sti^fyd  Uai.v.ersky 


IR/UV  Radiation  Conversion 
in  Optically  Pumped  CO-Ar  Plasma 

Relative  population 


Pump 

Laser 

IN 


V-V  exchange 


V-E  transfer  / 
UV  Radiation  OUT 


Vibrational  Quantum  Number 


Nonequilibrium  Thermodynamics  Laboratories 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


GRADIENT  DISTORTION  EXPERIMENT 

•  Transverse  temperature  gradient  is  created  by  absorption  of 
CO  laser  radiation  (1-2%  CO  is  added  to  the  flow  to  enable 
absorption  and  as  a  thermometric  element). 

At  P=100  Torr,  dT/dr  -  100  K/cm 

Routinely  used  in  CO  and  NO  optical  pumping  kinetic 
experiments. 

•  Plane  shock  wave  created  by  a  capacitor  bank  discharge. 
Measurements  oft 

•  flow  density  distribution 

•  Gas  temperature  and  transverse  temperature  gradient  before 
the  shock  arrival 


Radial  temperature  profile  is  parabolic: 
T(r)  =  Tc  -  (Tw-Tc)(r/R)2 


Nonequilibrium  Thermodynamics  Laboratories 


Non-Equilibrium  Tliermodynamics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 
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Non-Equilibrium  Tliennodynamics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 
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